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I ■ We study mass-radius relations for compact stars employing an equation of 

2 ■ state (EOS) of dense matter based on a Skyrme fluid. The zero-temperature 

^ ■ mean-field model is based on mesonic excitations, incorporates the scale break- 

er ■ 

ing of QCD, and accommodates baryons (nucleons) which arise as a solitonic 
configuration of mesonic fields. Stable configurations are obtained for central 
densities Pc/Pn < 5.0 where p„ = 2.575 x 10^^ g/cc is the nuclear saturation 
density. These Skyrmion Stars are mostly fluid, with a crust which we describe 
by the EOS of Baym, Pethick and Sutherland. Their masses and radii are in the 
range 0.4 < M/Mq < 3.6 and 13 km < i? < 23 km, respectively. The minimum 
spin period is computed to be between 0.7 ms < P < 2.1 ms. They appear 
to have a mass-radius curve quite different from either neutron or quark stars, 
and provide a suitable description of the heavier mass neutron stars discovered 
recently due to the inherently stiff EOS. Within the same model, we compute 
the dominant neutrino emissivity in neutron-rich /3-equilibrated matter, and de- 
termine the cooling behaviour of Skyrmion stars. 



Subject headings: skyrmions, compact stars, neutrinos 
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1. Introduction 

An oft-stated goal in the study of compact stars (neutron stars or possibly even quark 
stars) is the determination of the phase of matter in their interior through astronomical 
observations (Weber 1999, Lattimer 2000, Heiselberg & Pandharipande 2000). An equation 
of state is used to describe this phase; since it is an input into the equations that describe the 
structural properties of the star, it can therefore be constrained by stellar mass and radius 
observations (van Kerkwijk 2001, Lattimer & Prakash 2001, Lattimer & Prakash 2004). 
The challenge in this field arises, on the one hand, from our lack of theoretical knowledge 
about the EOS of matter at the highest densities, and on the other hand, from mass and 
particularly radius determination of compact stars (Haensel 2003, Heiselberg 2001). Despite 
improved many-body physics of nuclear and neutron matter of relevance to neutron star 
interiors (Schwenk 2004, Schwenk 2003), and increasingly accurate estimations of mass and 
radius (Manchester 2004), a number of EOS remain in use, with no definite or universal 
answer as to which is the best one at supranuclear densities. That being said, theoretical 
studies of the nuclear equation of state can, to an extent, be gauged for their relative success 
when confronted against data (particularly, accurate radius measurements) , and we refer the 
reader to the article by Lattimer & Prakash 2001 (and references therein) for details. 

Until recently, timing measurements of radio binary pulsars yielded accurately measured 
neutron star masses in the range 1.26 to 1.45 Mq (Thorsett and Chakrabarty 1999). Newer 
observations of binary systems of pulsars and white dwarfs allow for a larger range of masses; 
for example, the binary component PSR J0751-I-1807 with 2.11q;5M0 at 95% confidence 
(Nice et al, 2005). We also point out the possible confirmation (at 95% confidence) of at 
least one pulsar with M > I.68M0 by Ransom et al, (2005), following the recent detection 
of 21 millisecond pulsars in the globular cluster Terzan 5, 13 of which are in binaries. These 
observations of presumably heavier mass neutron stars seem to favor a stiffer equation of state 
than has been explored in the literature, at least for the central densities that we consider in 
this work (pc < S.Op^)- Intriguingly, it has been shown in Ouyed & Butler (OB98) that an 
equation of state based on the solitonic Skyrme model yields quite stiff equations of state and 
neutron star masses greater than 2Mq for standard radii of 15 km or so (while still obeying 
causality limits). In this letter, we extend the work done in OB98 in two main directions. We 
construct the equation of state for the Skyrme fiuid including isospin breaking effects in order 
to address asymmetric matter. This is more realistic than the approach in OB98, where only 
symmetric matter and pure neutron matter was considered, since neutrons begin to drip out 
of nuclei already at sub-nuclear densities, and the neutron star is generally characterized by 
an asymmetric mixture. We contrast the resulting mass-radius curve with that of neutron 
and quark stars constructed out of more conventional equations of state. We also address 
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neutrino cooling within the same model to obtain a more complete picture of the behaviour 
of Skyrmion stars, and compare our cooling curves to temperature data. Our main goal 
in this paper is to motivate, analyze and elucidate the consequences of using the Skyrme 
model as a description of compact star interiors, particularly in light of the recent discovery 
of neutron stars with masses greater than the canonical IAMq. This paper is organized 
as follows: In §2, we consider observational hints for large mass neutron stars. In §3, we 
describe the Skyrme model, its extension to dense matter at nuclear densities and beyond, 
and obtain an equation of state. In §4, we display the mass-radius curves that are obtained 
from this EOS, and consider dominant channels for neutrino emission and cooling of the 
star. Our findings and conclusions are detailed in §5. The Mathematica routines used in 
constructing the EOS for our model, as well as those used to generate cooling curves, can be 
accessed on the public domain http://www.capca.ucalgary.ca and downloaded for personal 
use. 



2. Observational motivation 

We are closing in on neutron stars both observationally (masses and radii) and theo- 
retically (EOS at high density). Observations of X-ray binaries offer the best support for 
larger masses with values of 1.8 - 2.4Mq (QuaintreU et al, 2003), 2.4 ± 0.27Mo (Clark ct al, 
2002) and 2.0 — 4.3M0 (Shahbaz et al, 2004) being reported. Although these measurements 
are not as precise as for radio pulsars ^ they seem to hint at the existence of heavy neutron 
stars. 

Such massive NSs are further suggested indirectly by (i) a possible interpretation, within 
current uncertainties in age and surface temperatures, of some colder neutron stars such as 
Vela and Geminga pulsars as moderately heavy stars (~ 1.5 — 1.8A'/q). These suggestive 
calculations are based on models of neutrino cooling taking proton and neutron superfluidity 
into account (Kaminkcr et al. 2001, Gusakov ct al. 2005). It should be noted, however, 
that an interpretation of these pulsars as low mass pulsars, being consistent within the same 
uncertainties, cannot be ruled out (Page et al. 2004). (ii) modeling of neutrino-driven winds 
from very massive (M ~ 2.OM0) and compact {R ~ 10 km) neutron stars which reproduced 
the r-process abundance pattern with reasonable success (Otsuki et al, 2000; Sumiyoshi 
et al, 2000), due to the short expansion timescale and high entropy that are required to 
obtain a high neutron-to-seed ratio. However, such a large mass can be criticized from the 



^We should note that there is only one NS in a binary radio pulsar system, the pulsar J0751+1807 for 
which a large mass 2.1 ± 0.5Mq has been derived (Nice et al. 2004). 
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preponderance of observational data on neutron star mass; (iii) recent studies by Popov & 
Prokhorov (2005) addressing the issue of existence and formation of massive Skyrmion stars. 
They considered different channels for the formation of massive rapidly rotating neutron 
stars using a population synthesis code to estimate numbers of massive neutron stars on 
different evolutionary stages, and found that a significant increase in neutron star mass due 
to accretion is possible for certain values of initial parameters of the binary. They concluded 
that, within the Skyrmion star scenario, rotating neutron stars with M > 2Mq can arise 
from binary evolution. The authors argued that a significant part of such heavy stars can 
be observed as millisecond radio pulsars, as X-ray sources in pair with white dwarfs, and 
as accreting neutron stars with very low magnetic fields. These massive Skyrmion star 
candidates (with the above mentioned observable signatures) can hopefully be tested in the 
near future; we refer the interested reader to Popov & Prokhorov (2005) for more details. 

The evidence for large neutron star masses is still being debated and one awaits confirmation. 
However, the large masses hinted at above pose a challenge to models of neutron stars based 
upon modern EOS even as uncertainties are being reduced by improved two- and three- 
body forces, relativistic effects and many-body calculations (Heiselberg &: Pandharipande 
2000). Even the stiffest EOS developed so far seems to be facing difficulty in accounting for 
the extreme values (up to 2.4Mq; sec discussion in Heiselberg 2001). This brings us to a 
theoretical discussion of the Skyrme model as an alternative. 



3. Theoretical motivation 

3.1. The Skyrme model 

Well before the discovery of Quantum Chromodynamics (QCD) as the fundamental 
theory of the strong force, Skyrme (1962) provided a model based on pionic excitations, 
in which the nucleon arises as a classical configuration of the pion fields with a non-trivial 
topology and a localized energy density (a soliton). This is remarkable in that it allows 
us to consider fermions (nucleons) in a bosonic (pionic) theory without introducing any 
exphcit fermion fields. The Skyrme model has been placed on a sound theoretical footing 
in subsequent analyses (Wittcn 1979) and it is now accepted as a valid description for 
certain aspects of the strong force. Its connection to phenomenological nuclear forces is 
well established (Nyman and Riska 1990), and its justification based on large arguments 
('t Hooft 1974) stemming from a quark-gluon picture indicates that it acts as a "bridge" 
between QCD and nucleon-nucleon (n — n) forces. Several detailed reviews on the Skyrme 
model may be found in the literature (Holzwarth and Schwesinger 1986, Zahed and Brown 
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1986, Schechter and Weigel 2000). 

3.2. The Skyrme model and gravity 

Skyrmions coupled to gravity have been previously explored in Glendenning et al. 1988, 
Droz et al. 1991, Bizon and Chmaj 1992, with the conclusion that Skyrmions with large 
topological charge N cannot be models of neutron stars due to gravitational collapse at micro- 
sizes {N ~ 10-*^^) compared to neutron stars {N ~ 10^^). Furthermore, they are unstable to 
N — 1 configurations if non-spherical shapes are disallowed. Therefore, we adopt a more 
conservative approach in which the star is composed of an interacting ensemble (at the 
mean-field level) of iV = 1 skyrmions (unit baryon number). In the spirit of nuclear mean- 
field approaches to the equation of state, we utilize the skyrme model for the microscopies, 
and couple the ensemble of = 1 skyrmions to gravity through the stress-energy tensor 
T^^ of the dilute skyrme fluid (no overlap between skyrmions). This implies that gravity 
becomes strong only at very large N (unlike in afore-mentioned studies where the entire star 
is treated as a topological object coupled to gravity) and that instability towards single- 
particle emission is not relevant in our model. The skyrmion profile and its mass-radius 
characteristics are very different than obtained in previous works. This approach will provide 
for a much more realistic description of neutron stars, while maintaining a connection to the 
Skyrme model. 

3.3. The Skyrme fluid 

For astrophysical applications, one needs a many-body description of the nuclear force, 
and for neutron star densities, an in-medium approach based on the relevant degrees of 
freedom in the theory. As a successful and theoretically sound model, the Skyrme Lagrangian 
bears further investigation as regards its application to astrophysics. A particularly useful 
mean-field model was developed in the work of Kalbermann (1997), who studied Skyrmions 
in a dilute fluid approximation, where auxiliary dilaton (cr) and vector (cu) flelds couple to the 
Skyrmion and carry the information of density and temperature. The (non-zero) expectation 
values of these background flelds are determined from the Skyrme Lagrangian itself, in the 
spirit of the mean held approach. 

The Lagrangian for the Skyrme model, augmented by the a, u fields, and including isospin- 
breaking effects from the p meson as well as explicit scale-breaking effects from the dilaton 
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and quark masses is given by 

F2 1 

rr iTrrr.„r/^ + 

32e' 



e^^ 



P2 2 

+ e2'^^-^Tr(C/ - 1) - ^vc^^S'' - B{\ + e^"(4(7 - 1)] 



1,. . .0 1 o „ 1 

4^ 



where 



L^^U^d^U . (2) 



In eqn.(l), F^^ is the pion decay constant, Cg is the Skyrme parameter, and g\ is the uj N 
couphng strength. The a;-meson couples to the topological (baryon) current 

B^=^TAiU^d^U){U^d^U){U^d,U)] . (3) 
The p meson, introduced here to extend our model to asymmetric matter, couples to the 

— * 

isospin current I'^ with coupling Qp, and its field strength is non-abelian F^^ — df^pi, — di,Pfj_ + 
9p{Pti ^ Pu)- The p meson is viewed here as a massive Yang-Mills field, and the model is 
rcnormalizable with the introduction of a Higgs field. To obtain the equation of state, wc will 
work in the limit of infinite Higgs mass, which simplifies the equations of motion for the mean 
fields (Scrot and Walccka 1986). This is a standard approach in p, uj, cr-typc Walccka models 
(Quantum hadrodynamics) which closely resembles our model. Scale invariance is broken 
by quantum fluctuations of the QCD vacuum and flnite quark masses, and this feature is 
manifest since the dilaton mocks up the scale anomaly of QCD (Nielsen 1977, CoUins et al, 
1977) 

g 

where is the trace of the energy-momentum tensor, D^^ = T^^x^ is the dilatation current, 
as the strong coupling constant, and G^j, the gluon field strength. The sum over q runs over 
the up and down quark fiavors, and we have ignored the anomalous dimension. Consistency 
with the trace anomaly demands that the quark condensate behaves as {qq) {qq)e^'^ under 
a scale transformation r — > re~^. Its connection to eqn.(l) is through the Cell-Mann Oakes 
Renner relation 



F^ml = -2(m„ + ma){qq) . (5) 
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Prom the above Lagrangian, it is easily verified that the gluonic contribution 

T; = (roe'^)^ ■ (6) 

In more famihar terms, Tg/e — —4B, where B is the bag constant. The dilaton field also 
provides the required intermediate-range attraction in the A^A^ force. The a;-meson provides 
the right sign for the isospin-independent spin-orbit force. It can also stabilize the soliton in 
the absence of the Derrick term (£4) (Adkins and Nappi 1984a). For the sake of completeness, 
we have included the heretofore ignored pion mass term, which has recently shown to be 
important for multi-skyrmion systems (Battye and Sutcliffe 2005). In fact, these authors 
find that alternate minima of the Skyrme energy functional can arise if either the baryon 
number or the pion mass is large. For the (physical) value of the pion mass, and the product 
ansatz for the Skyrmions that we use in this work, we can continue to use the "hedgehog" 
ansatz ^ 

U = exp(if.fF(r)) , (7) 
with T the Pauli matrices. With some algebra, one can show that the Skyrmion mass is 



M 



Moe'^° + M^N] Mo = 4n j dr r^Mo{r) ; 
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+ 



sin^F 



nmlF^ J r^(l - cosF) , 



sin^F .2 



(8) 
(9) 



where (Tq = (a), ujq = {uj) and po = (p) henceforth denote mean-field values. The 
resemblance of the above form to the free Skyrmion case comes from the association of a 
with broken scale invariance. With the product ansatz for the Skyrmion 



UB=N{r, Ri, .., Rn) = U{r - Ri)..C/(r - Rn) , 



(10) 



and the sum ansatz for the cr,u!, p fields 



crB=N ^ CTi + .. + aN] ujn ^u;i + ..+u;n; Pn^Pi + -- + PN: 



(11) 



^The product ansatz is used for simplicity. For truly multi-skyrmion shell-like solutions with large baryon 
numbers of 20 or more, the pion mass plays an important role. 
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the mean-field treatment implies that the equations for the Skyrme fluid may be derived 
by solving the Eulcr- Lagrange equations from L for a single free Skyrmion in the absence of 
background fields, and rescaling as 

r^e^^V, a;^e^°a;, p ^ e^V (12) 

in the presence of the medium (Kalbermann 1997). The free Skyrmion mass Mq and the 
contribution from explicit scale breaking (the pion-nucleon sigma term M^at) can be formally 
obtained by noting that the solution for F(r) is (Adkins and Nappi 1984b) 

Fir) = 0.759 . (13) 

r 

The "effective mass" of the Skyrmion in the fluid is then M. We now have the tools to 
construct the equation of state for the Skyrme fluid. 



3.4. Equation of State 

To obtain the equation of state, we require a relation between the pressure and energy 
density of the fluid. In the mean field approximation, the energy of N Skyrmions per unit 
volume is given by 




- -e^''°m^u;o + Qv^^pv 

1 1 

- -e^'^^mjp^ + -5pPoP/ 

(14) 

where'^ Py is the net (baryon minus antibaryon) density and pi is the net (proton minus 
neutron) isospin density. We choose nip — 770 MeV, and we will fit the binding energy for 
symmetric matter (-16 MeV per nucleon) and symmetry energy (32 MeV) at saturation to 
obtain the parameters gy and gp respectively. Here, the particle and anti-particle distribution 
functions are 



"^the extra 1/2 factor in the p coupling compared to a; coupling is due to the conventional definition of 
the Pauli isospin matrices 
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= (exp^ + ' ; = (^exp^ + ' (15) 

with the dispersion relations ef°* = EP'^^ + gvoUo + gppo/'^, e^^"* = ^p^^Hffyc^p -gpPo/2, 
^prot ^ £;prot_^^^^ + ^^^^/2, e^^"* = E^^^'-gvUJo-gpPo/2 where £;p™* = = ^/^ + JP 
(our model yields the same Dirac effective mass for neutrons and protons). The chemical 
potentials satisfy 



^liPFi) ^ IJ'i; prot, neut (16) 

where pp^ is the Fermi momentum of species i. The quantity Ep is related to the energy of 
an individual Skyrmion Est by a Lorcntz boost of the static solution (Kalbermann 1997). It 
implicitly includes a contribution from the pion mass term from Eq.(9). 

Since the scale of thermal excitations in the Skyrme fluid is much less than the baryon 
chemical potential, we will work at T = for this section. ^ The Fermi distributions reduce 
to theta functions, and we find 



\ ^ iio^uy \i x-prot ^neut' / 1 ^\ 

3^^ ' 3^^ ■ ^''^ 

The equations of motion for the mean fields are obtained by minimization of the energy 
eqn.(14) with respect to the a, cu and p fields : 



dEv 2 / /-Pw \ fdEp\ dV{ao 



dao (27r)^ \J J J \dao J dao 

- e^''°mlujl-e'"''mlpl . (18) 
dE 

= —^^-e'"''mluj^ + gvpv ■ (19) 
dujQ 



= ^ = + ?f P. . (20) 



The symmetry energy can be calculated in our model to be 



"^While considering neutrino cooling, we naturally admit finite temperatures. 
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where pp — (37r^py/2)^/^. At saturation py = — 0.154 fm~^ and with an effective mass 
oi M — 600 MeV and a symmetry energy of 32 MeV ^, we get Qp — 7.70. In asymmetric 
matter, beta-equihbrium demands Pneut — A*prot + A*eiec, while local charge neutrality dictates 
that riprot = neiec- These conditions can be substituted back in the mean field equations as 
follows : 
Let 

X = P^; = -^, (22) 

PFprot ^nevLt + nprot X^^ + 1 

The energy Ey can be expanded about the symmetric point y = 1/2 to a good approximation 
as 



Ev{py,y)^ Ey{pv, 1/2) + as{py,ao{py)){l-2yf + ... , (23) 

where ... arc ignorablc corrections. Studies of neutron matter show that this expansion 
works well over a wide range of values for y and py (Prakash et al. 1988). Then, beta 
equilibrium along with charge neutrality implies 

(37rVz/)'/' = 4a,(pv, ao(py))(l - 2y) , (24) 
which fixes a solution y = yo for any baryon density. It follows that 

PFp.„t = (3vrV2/o)'/'; PF_t = (3vrV(l-2/o))'/'; pi = pv{2yo - I) ■ (25) 

Performing the momentum integrations in eqn (18) analytically, the equations of motion 
become 



^The symmetry energy in dense matter determines the threshold for the direct urea process involving 
protons {n p + e~ +De)- This is an extremely efficient neutrino cooling mechanism, and will be addressed 
in the section on cooling. 
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= 



+ 

= 
= 



M2 
M2 



M2 



dcTQ 

dEy 

dujo 

dEy 



-e'^'^°mluJo + gvpv ■ 



(26) 
(27) 

(28) 



For symmetric matter, a fit to the binding energy By = Ey/pv — M of -16 MeV at saturation 
yields gy = 7.54 (OB98). In our case, for neutron-rich matter (|/ -C 1), this value leads to 
a positive By (see Figure 1), implying that the matter is unbound, as expected. We do 
not expect gy to be sensitive to the isospin composition, as the a;-meson couples to the net 
baryon density. Adopting this value of gy, the above equations can be solved for the mean 
fields to yield the energy per baryon from eqn (14) at any baryon density. 

The pressure of this ensemble is given by 

opv 

The respective equations of motion show that the contribution of the vector meson u to the 
pressure grows with density and is positive, a consequence of acting in the repulsive channel. 
The dilaton a gives a negative pressure and tends to decrease the (magnitude of) binding 
energy. The p meson acts to stiffen the equation of state, particularly at higher densities. 
These features drive the pressure-density and mass-radius relations that we obtain in the 
following section. 



3.5. Fit pcirameters for neutron-rich matter 



In order to reproduce neutron matter phenomenology, four free parameters (oj, i = 1, 4) 
are added to the conventional dilaton potential.^ Then, 



B [l + e^'^(4cT-l)]+5[ai(e-'^-l)+a2(e'^-l) 

+ a^ie^"" - 1) + a4(e=^" - 1)] , (30) 

with the "bag constant" B ~ (240 MeV)^. Terms of the form e""' are the only ones that can 
be added to the potential in accordance with the anomaly condition 

at cr = 0, implying a\ = 2a3 + 304/ The remaining three parameters of the dilaton 
potential are constrained by demanding that at saturation density (p„ = 0.154 Baryons/fm^), 
the effective mass M/Mq = 0.60, the compressibility K — 270 MeV, and that eqn.(26) is 
satisfied. On fitting to these properties, we find 



ai = -29689.7, = -274771.0, (32) 
ag = 222839.0, 04 = -66865.6 . (33) 

Although large in comparison to the fitting parameters obtained in OB98, there are 
large cancellations that lead to typical values V(a^ ~ 0.06-B for the entire range of a (and 
hence density) explored in this work. Figure 1 shows the binding energy per nucleon as a 
function of the density as well as the pressure-density relation (EOS) and the sound speed. A 
comparison to more conventional equations of state, as well as the EOS of OB98 is presented 
therein. The increased stiffness with respect to the results of OB98 is evident. It is due to the 
limited attraction provided by the dilaton, and the extra repulsive effects from the p-meson. 
We now turn to apphcations of the Skyrme fluid EOS. In the next section, we discuss the 
global structure of Skyrmion stars and their coohng properties on the same footing. 



*'Thcre is also the possibility of modifying the w field potential (see Kalbermann (1997) and references 
therein) . 

^One chooses the term multipHed by a\ to have a negative power of a in order to avoid the introduction 
of a second minimum in the potential for cr < 0. The only sensible minimum then remains the one at cr = 0. 
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4. Implications of the Skyrme EOS 

4.1. Mass-radius relations 

We proceed now to construct models for stars using the EOS developed above. This 
is done by integrating the general-relativistic equation of hydrostatic balance, including 
rotation, upto the surface of the star where the pressure vanishes (Tolman 1934; Oppenheimer 
& Volkoff 1939). For the rotating configurations we use the computer code RNS (Stergioulas 
& Friedman 1995; 1998). We specify the equation of state and the central energy density, 
and the code computes models with increasing angular velocity until the star is spinning with 
the same angular velocity as a particle orbiting the star at its equator (see Ouyed 2002). 

Figures 2a and 2b show the resulting stellar masses and radii as a function of the central 
density, respectively. For the rotating case, the masses are in the range 0.4 < M/Mq < 3.6 
while the radii (equatorial circumferential radius; [proper equatorial circumference] /27r) are 
calculated to be 18.6 km < R < 23.0 km. In Figure 2c, we show the resulting Mass-Radius 
plane. Note that the maximum mass configuration does not correspond to the maximum 
radius case. In general, for a given central density, rotation allows the masses and radii to 
increase by 30 % and 40%, respectively, when compared with the non-rotating cases. The 
amount of mass in the outer region of the star (the crust region - p < p„ - constructed 
using the EOS of Baym, Pethick & Sutherland, 1971) decreases with rotation. The crust of 
rotating Skyrmion stars constitutes less than 5% of the total mass while it averages 20% for 
the static configurations (ref. Figure 5 of OB98). The minimum spin period for Skyrmion 
stars using our improved EOS was computed to be 0.7 ms < P < 2.1 ms not very different 
from what was found in Ouyed (2002). 

A comparison of the mass-radius curve for the Skyrmion star with that of neutron and quark 
stars (see Figure 3) shows a qualitative difference that arises from the physics of the dilaton. 
The limited attractive force provided by the dilaton, and the inclusion of the p meson leads to 
large stiffness even at moderate central pressures when compared to ordinary neutron stars. 
The radius of Skyrmion stars increases with mass within a certain mass window, despite 
increased gravitational forces, similar to the curve for self-bound quark stars (see rescaled 
curve in Fig. 3) and neutron stars described by a highly stiff equation of state (see curve 
labelled MS) based on a relativistic mean-field approach (eg. Mueller and Serot 1996). 
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4.2. Neutrino cooling of the Skyrme fluid 



A unique method to probe the internal composition of a neutron star is by tracing 
its temperature evohition with coohng simulations (Pethick 1992, Page 1998, Yakovlev et 
al. 2003). Neutron stars are born in supernova explosions and as the neutron star cools, 
neutrinos begin to free stream and essentially leave the star without further interaction. As a 
consequence, about 30 seconds after the explosion, the cooling of neutron stars is controlled 
by neutrino emission. Importantly, the neutron stars remain luminous enough during the 
neutrino cooling era, so that the surface temperature can be extracted from space telescope 
data and the theoretical cooling curves can be confronted with observations. There is a host 
of well-known neutrino emission processes that operate in the crust and core of the neutron 
star (Yakovlev et al. 2001), and which one dominates the cooling depends mainly on the 
temperature and to a lesser extent, on the density. Since cooling rates are dependent on 
neutrino emissivities, which in turn depend on the dense matter equation of state employed, 
it is instructive to compute the dominant neutrino emission rates. We now proceed to do so 
for the Skyrme fluid. We address only the global thermal energy balance of the star during 
the neutrino and photon cooling eras, so that we may employ 



to compute the cooling history. The required inputs to the above equation are the total 
specific heat at constant volume Cy, the neutrino emissivity which determines the neutrino 
luminosity as ~ ^uV , where V is the cooling volume of the star and L-y, the photon 
luminosity which determines the late time (> 10^ — 10^ yrs.) cooling behaviour. We will 
analyze this equation by approximating — 0, which defines the neutrino cooling epoch, 
or L,^ = which defines the photon cooling epoch. In asymmetric matter, the direct urea 
process (in medium /3-decay) n — > p + e" -|-Pe and p + e — > z/g as well as the modified urea 
process n+n n+p+e+z/g and n+p n+n+e+u^ become important channels of neutrino 
emission. The former is kinematically allowed only if energy and momentum conservation, 
expressed by the triangle inequality Q{pFe +PFp —PFn) is satisfied. This happens at a critical 
value of y. On using the above condition, the threshold where direct urea can start to occur 
is (Lattimer et al. 1991) 




7 



(34) 



Z/crit — 




(35) 
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Since charge neutrality demands pp^ — PFp, this imphes that 



ycrit = g ■ (36) 



Only if y exceeds 1/9 can we have a direct urea process. The associated neutrino rate 
is (Lattimer et al. 1991) 



. 1/3 

dUrca 



Here, the neutron and proton Landau effective masses ^ are a function of density and 
are given by 
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Pk 



(38) 
(39) 



where pp^ = {'in'^ pyyY^^ ; pp^ = (37r^py(l —y)Y^'^. From Figure 4 for the proton frac- 
tion y{pv), it appears that neutrino emission by the direct urea process sets in for densities 
close to 2.5pn- The modified urea, in which energy-momentum conservation is facilitated by 
a spectator nucleon, occurs at all densities, and the associated neutrino rate is (Priman & 
Maxwell (1979)) 



^MUrca ^ 2.45 X lO^^ 

/ 



m: 



Q!Mi: 



GeV J \ GeV 
1 



, 1/3 
PV \ ,.l/3j.l 



Pn 



y^' "I g erg cm ^s ^ ; 



VI + 0.173 (^)'" 5^ 



(40) 



The calculation is described in detail in Friman & Maxwell (1979), who worked in neutron 
matter with a phenomenological description of nuclear forces, and explicit nucleon fields. In 



^The Landau effective masses m* and m* account for momentum-dependent interactions at the Fermi 
surface. 
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our case, we used the Skyrme model to derive eqn.(40). The Skyrme model can reproduce 
the long-range part of the A^A^ interaction, namely, the one-pion exchange (Nyman & Riska 
1990), which is the most relevant piece for the modified urea process at these densities. 
In obtaining the above result, eqn.(40), the following points should be noted. We used 
the pion-nucleon coupling constant value Qt^nn = 11.9. This is the predicted value in the 
model (Adkins and Nappi 1984b) that we are using, compared to the actual experimental 
value of 13.2. The model parameters are chosen to reproduce the nucleon (938.9 MeV) 
and Delta (1232 MeV) mass. Although we have used a relativistic model for the nucleons 
thus far, we approximated their dispersion relations by their non-relativistic counterparts for 
calculational simplicity. Thus, Ep = p^/2M — p"^ /^M^-\- .. . The relative error in retaining 
only the leading term is about p|n/4M^ ~ 10% if we consider densities as high as 2p„. 
This does not affect our final result for the emissivity in any considerable way. At higher 
densities, the Fermi momentum of the nucleon is not small compared to its effective mass M, 
and a fully relativistic treatment is required. Here, we restrict our application to Skyrmion 
stars with central densities pc < 2p„ so that our non-relativistic approximations hold. This 
restriction to low density also means that the direct urea process remains forbidden. Prom 
the plots for mass versus central density and mass versus radius, we observe that this limits 
us to M < 2Mq and R < 15km. Thus, we are well within current observational hmits. 

The neutrino luminosity Lj, is now determined as 



where R denotes the radius of the star. The specific heat of the baryons can be expressed as 



where /(ep) = (1 + exp(^^^))~^ and we have omitted a term that is negligible in the 
degenerate limit. In this limit, one can extract the density of states D[ep) at the Fermi 
surface and express the specific heat as 




(41) 




(42) 




(43) 



The total specific heat is then 




(44) 
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With eqns.(41) and (44), we can now proceed to compute the cooling curve for the star. In 
fact, an analytical estimate is easily made (Page et al. 2004). Performing the integrations 
in eqns.(41) and (44), we obtain 

= NT^] Cv^CT . (45) 

where C and N depend on density but not on temperature. The cooling equation (34) can 
then be integrated to yield 

T^(^)"\-'/« (46) 

for long-term cooling in the neutrino dominated epoch. The "internal" temperature 
T is, strictly speaking, the temperature at the bottom of the crust, and can be related to 
the temperature at the surface of the crust via the so-called Tsuruta law (Tsuruta 1979) 
which is an interpolation between the surface and the crust temperatures. The formula 
for thick crusts, as is the case in our model, is (see also Shapiro & Teukolsky 1983, p330) 
Tb = {lOT^f/^; the subscript "m" stands for mantle so that T = T^, the subscript "b" 
for the bottom of the stellar envelope (which is also the top of the crust). It is likely, 
however, that very soon after the neutron star's formation, the crust becomes supcrfiuid, 
and is therefore almost isothermal. We can then obtain the surface temperature Tg at the 
top of the stellar envelope (atmosphere) directly by employing the Tg — Tf, relation, first 
derived in Gudmundsson et al. (1982), and utilized for practical cooling simulations in Page 
et al. (2004). With n = T, we find 

T, = 0.87xlO%g,J'/\Tsf'+'' . (47) 

Here, g^s^^ is the surface gravity gs measured in 10^^ cgs units (see eqn.(50)), Tg is temperature 
in units of lO^K and a parameterizes the composition of the envelope (we choose a — 0.05 
which is typical of an iron envelope). 

The effective temperature as seen by a distant observer is then = Tpy^l — Rsch./R where 
-Rsch. is the star's Schwarschild radius. The resulting cooling curves {T^ versus time) are 
plotted in Figure 5 for Skyrmion stars with central density pc = 2p„. The cooling curves 
for the non-rotating configurations correspond to M = 2.2Mq, R — 1& km and those for 
the maximally rotating ones correspond to M = 3.OM0, R — 22 km. Note that we included 
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photon cooling which starts to dominate when Lj ^ L^,. In that era, the T vs. t behavior is 



C 



l/4a 



1 



T=(^) r-, (48) 



where C and a are the same constants as before, and 



S - 4.08 X 10i^-32a^2^ g^^^ ergs/sec K^+^" ; (49) 
M e^" 

g,,, = 1.334—-^; (50) 
1 



0-296^,^ 



(51) 



with RiQ the stellar radius in units of 10 km. The photon luminosity in units of ergs/sec 
is given by = 5'T^+^°. Irrespective of the initial temperature {T(t = 0)), we observe 
that Tg°° in Figure 5 traces the temperature data, which arc mostly constrained only by 
upper limits. Photon cooling starts to dominate close to a million years after the neutron 
star's formation. The exceedingly rapid cooling of the star then makes it thermally invisible 
(Tg°° below lO^K). The parameter accounts for a possible suppression of the nucleonic 
specific heat by superfluidity. The electrons then control the specific heat, given by Cy — 
IQ-o^s (.^pi^^sapped) ^.^j^ ~ 2 a typical value (Page et al., 2004). We have shown curves 
for two cases, ag — 1.0, 2.0. Note also that for the rotating case, neutrino cooling is more 
efficient but photon cooling takes over at a later stage as compared to the non-rotating case. 
This is due to the smaller crustal mass fraction for the rotating case. Our estimates of the 
cooling rate are bound to be affected mildly by inclusion of other neutrino emission processes 
such as neutrino bremsstrahlung and Cooper pair-breaking and formation processes, as well 
as effects from time-evolving stellar envelopes, but our results capture the essential cooling 
history of a Skyrmion star. 



5. Conclusions 

In this work, we have examined the case for considering the Skyrme model as a viable 
description of dense nuclear and neutron-rich matter, with direct application to neutron 
star interiors. There are clear indications that conventional neutron matter equations of 
state cannot explain the observation of some very heavy neutron stars, unless the central 
densities are assumed to be fairly high. The Skyrme Lagrangian, within the mean-field 
approach employed here, and with the inclusion of the dilaton and vector mesons, serves to 
generate quite stiff equations of state that can describe the structure of these heavy stars. 
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Importantly, the Skyrme model incorporates features of QCD such as scale breaking (the 
trace anomaly) and chiral fields that map to low-lying excitations of the vacuum. It should 
be noted, however, that the use of the dilaton as an interpolating field to mimic the (large) 
scale-breaking in QCD appears at odds with the estimated glueball mass of 1.5 GeV and 
its stiffness in matter (Birse (1994)). A better in-medium description of the Skyrme model 
would have to take this fact into consideration. 

The mass-radius curves for Skyrmion stars are distinct from ordinary neutron stars, essen- 
tially due to the qualitatively different nature of the dilaton potential, which respects the 
scale anomaly, and is fit to nuclear matter properties. We have described both rotating 
and non-rotating configurations for Skyrmion stars, the essential difference between them 
being the crustal mass fraction, and the total mass that they can support (rotating con- 
figurations can support considerably more mass). Their masses and radii lie in the range 
0.4 < M/Mq < 3.6 and 13 km < i? < 23 km, respectively. The minimum spin period lies 
between 0.7 ms < P < 2.1 ms. Skyrmion stars offer a possible explanation for observations 
of large mass neutron stars. Such stars are also indirectly hinted at in astrophysical models 
of the r-process, cooling curves of neutron stars, and have been examined in population 
synthesis studies. 

We have also studied cooling of Skyrmion stars by identifying the dominant neutrino emission 
process, the modified urea, and computing the corresponding cmissivity. The undcrpredic- 
tion of qa and g-n-NN in the model implies that the neutrino emissivity is also underestimated 
by a factor~ 5. However, this has only a shght effect on the cooling rate since the temper- 
ature depends on the l/6th power of the ratio C/N (see eqn.(46)). Along with the specific 
heat of degenerate matter, this allows us to sketch the cooling history of the star from tens 
of years up to a million years when the star rapidly cools to unobservably low temperatures. 
A more accurate cooling curve can be produced by taking into account the thermal conduc- 
tivity of the material, density dependence of superfluid gaps, and other neutrino emission 
processes. We have included photon cooling, assuming a generic envelope at the surface. 
While these surface structures are not direct consequences of our model, they are "pinned 
on" here for a more realistic picture, and to enable qualitative comparisons to the neutron 
star temperature data. 

An additional contribution to the neutrino fiux that we have ignored in this model comes 
from the elcctroweak decay of charged pions to neutrinos. Furthermore, neutral pions can 
also decay into ncutrino-antineutrino pairs. Although this process is forbidden by helicity 
conservation in vacuum, the presence of a medium implies that decay can occur in a frame 
moving with respect to the preferred rest frame. In fact, this process has been addressed 
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within the Skyrme model very recently (Kalloniatis et al. 2005). The reason we have 
neglected such "picnic" contributions is that the corresponding neutrino rate is penalized by 
an exponential factor of exp(— m^/T). Since T <^ 711^^ in our case, the neutrino emissivity is 
exponentially small. The dominant contribution in our model is indeed from the modified 
urea process, and from the direct urea process at much higher densities. 

Assuming a mean-field approach to the Skyrme fluid and ignoring Skyrmion overlap at 
high density is no doubt a simplification, but our aim in this work has been to gauge the 
model's main features, and examine its potential for the astrophysics of neutron stars. While 
improvements can surely be made on the theoretical aspects of the assumed model, we find 
that it displays qualitative features regarding its mass-radius curve that are very different 
from ordinary neutron stars. This could bear importantly upon the observations of heavier 
mass and/or large radius neutron stars as well as their cooling history. 
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Po/P„ 



Fig. 1. — a). Binding energy per nucleon in neutron-rich matter, as a function of the central 
density, b). Pressure versus density for our model (JO) compared to OB 98, BBB97 (Baldo 
et al. 1997), WW89 (Weber and Weigel 1989) and G89 (Glendenning 1989). c). Speed of 
sound ( in units of c) for our model as a function of density. 
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Fig. 2. — a). Gravitational Mass versus central density for zero temperature Skyrmion stars 
in hydrostatic equilibrium. Both rotating and static configurations for our model and OB98 
are shown to illustrate the differences, b) . Radius versus central density for Skyrmion stars 
in hydrostatic equilibrium, c). The Mass— Radius plane. 
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Fig. 3. — Mass-Radius relationship characteristic of neutron stars (NSs; e.g. Glendenning 
1997 (soft), Mueller and Serot 1996 (stiff)) and quark stars (SSs; e.g. Li et al. 1999) as com- 
pared to Skyrmion stars (SkySs). The SkySs model has been re-scaled for better comparison. 
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Fig. 4. — Proton fraction y versus baryon density in our model. Direct urea neutrino emission 
becomes allowed at y ^ 0.11 or p ^ '^.^Pn- 
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Fig. 5. — Cooling curve for Skyrmion stars with central density pc = 2pn, for maximally 
rotating (dashed) and static (solid) configurations, is a parameter related to crustal su- 
perfluidity (see text for details). Temperature data points are from Schaab et al. 1999. 



